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Two-photon loss mechanisms often accompany a Kerr nonlinearity. The kinetic inductance ex- 
hibited by superconducting transmission lines provides an example of a Kerr-like nonlinearity that 
is accompanied by a nonlinear resistance of the two-photon absorptive type. Such nonlinear dis- 
sipation can degrade the performance of amplifiers and mixers employing a Kerr-like nonlinearity 
as the gain or mixing medium. As an aid for parametric amplifier design, we provide a quantum 
analysis of a cavity parametric amplifier employing a Kerr nonlinearity that is accompanied by a 
two-photon absorptive loss. Because of their usefulness in diagnostics, we obtain expressions for the 
pump amplitude within the cavity, the reflection coefficient for the pump amplitude reflected off of 
the cavity, the parametric gain, and the intermodulation gain. Expressions by which of the degree 
of squeezing can be computed are also presented. 

PACS numbers: 42.50.Gy, 42.65.Yj, 42.50.Dv 
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I. INTRODUCTION 

Sensitive superconducting microwave devices such as 
SIS mixers Q, @ and parametric amplifiers @, have 
been devised which achieve performances close to the 
quantum limit. Phase-sensitive Josephson-junction para- 
metric amplifiers have been constructed whose noise per- 
formance exceeds that of the quantum limits imposed on 
linear phase-insensitive parametric amplifiers |5|. These 
phase-sensitive amplifiers have been used to generate 
quantum mechanical states of the electromagnetic field, 
called squeezed states, whose noise in one amplitude com- 
ponent is reduced below that of vacuum fluctuations. 
The kinetic inductance of superconducting transmission 
lines could also be used to make low-noise parametric am- 
plifiers. However, associated with the kinetic inductance 
is a nonlinear resistance that can degrade device per- 
formance. These nonlinear effects are relatively strong 
in superconducting striplines and microstrips due to the 
nonuniform distribution of the microwave current along 
the cross section of the transmission line. Along the 
edges, where the current density obtains its peak value, 
the current density can become overcritical even with 
relatively moderate power levels. As a result, the super- 
conducting current density may vary and, consequently, 
both inductance L and resistance R per unit length be- 
come current dependent according to the form |g 



L = L + AL 



(1) 



inductance provides a Kerr-like nonlinearity suitable for 
the construction of parametric amplifiers which employ 
four-wave mixing. The nonlinear resistance, to lowest 
order, is of the two-photon absorptive type. To aid in 
the design of parametric microwave amplifiers which em- 
ploy kinetic inductance we have preformed an analysis of 
cavity parametric amplifiers employing a Kerr nonlinear 
element for gain and a two-photon absorptive loss. Al- 
though the analysis was carried out with a specific appli- 
cation in mind Q, it is more generally applicable, since 
two-photon absorptive processes often accompany Kerr 
nonlinearities. There are optical [1 H E3 and mechani- 
cal j 1 1 1] systems with such combinations of nonlinearities. 

Squeezing in a parametric amplifier with a two-photon 
absorber has been studied by a number of workers 
[toUT^. ITiil Il4| . In the analysis provided here, we present 
expressions for the amplitude of the pump field within the 
cavity, the reflection coefficient for the pump off the cav- 
ity, the intermodulation gain, and the degree of squeez- 
ing. The first, second, and third of these quantities are 
particularly useful for extracting model parameters from 
experimental data. The equations of motion are derived 
using the input-output theory of Gardiner and Collett 
|l5l Il6| . The undepleted pump approximation is then 
made, allowing the pump field inside the cavity and the 
pump field reflected from the cavity to be calculated. 
The small signal response is then obtained by lineariza- 
tion about the pump field. 



II. THE HAMILTONIAN 



R = Ro + AR 



(2) 



where I (I c ) is the total (critical) current. The kinetic 



A lossless transmission line resonator having nonlinear 
kinetic inductance is discussed in appendix A and the 
effect of nonlinear losses associated with the kinetic in- 
ductance is discussed in appendix B. Here we consider 
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test port 



nonlinear resonator 

H = ft WqA^A 

+ (ft/2) KA^A^AA 



linear dissipation port 



a 2 

nonlinear dissipation port 



Y 3 f 



and 



H a 3 



(kv?Uxia\{uj)az{uj) 



(7) 



The linear coupling of the bath modes a\ and 0,2 to the 
cavity mode A is modeled by the hopping Hamiltonians 



FIG. 1: The model includes nonlinear resonator coupled to 
three ports, a test port, a linear dissipation port and a non- 
linear one. 



the case where the external signals employed for exter- 
nally driving the resonator are all close in frequency to 
one of the resonances at lu . As we discuss in the ap- 
pendix, under some conditions, which are assumed to be 
satisfied, all other modes of the resonator can be disre- 
garded. In this case the Hamiltonian of the nonlinear 
resonator can be written as 0, 0] 



H r = Huj A^A+ -KA^A^AA, 



(3) 



and 



Hn 



(8) 



(9) 



The two-photon absorptive coupling of the resonator 
mode to the bath modes 0,3 is modeled by a hopping 
Hamiltonian in which two cav ity pho tons are destroyed 
for every bath photon created [20. 12H I2H I2I 



H Ta =H dcu^A^ A^ a 3 {u) + K&al(u})AA] 



(10) 



All the modes in this model satisfy the usual boson com- 
mutation relations. 



where the Kerr constant K is given in Eq. IA36I 
As seen in Fig. 1, the resonator is coupled to a test port 
(labeled as a\) serving as the input-output port. Oper- 
ated as an amplifier, the signal returned or "reflected" 
from the input port is larger than the incoming signal. 
This mode of operation, at microwave frequencies, is re- 
ferred to as the negative-resistance reflection mode. Two 
extra fictitious ports are added in order to theoretically 
model dissipation [T3. Port 0,2 serves as a linear loss 
port. Port a 3 serves as the two-photon loss port. The 
coupling of the a 3 loss mode to the resonator mode A is 
nonlinear and given in Eq. 

It is convenient to write the Hamiltonian as a sum of 
terms, 



H — H r + H„ 



H n 



H 



H 7 



Ht 3 , (4) 



each representing the Hamiltonian for a component of 
the system. 

The three ports coupled to the resonator (see Fig. 1) 
serve as baths. One bath models the external modes 
that couple to the resonator mode through the port that 
serves both as the input port and as the output port. 
The Hamiltonian H a \ for this bath is given by 



H a i = / du}hwa\{uj)a\{uj) 



(5) 



The other two baths are associated with the linear and 
nonlinear cavity losses and their Hamiltonians are given 
by 



H a 2 = J duhuja^uj) 122(0;) 



(6) 



III. THE EQUATIONS OF MOTION 

Since the creation and annihilation operators appear- 
ing in Eqs. (??) through (|10|) do not have an explicit time 
dependence, the Heisenberg equation of motion for these 
operators has the form 



«§-p,m 



("I 



where H is the total Hamiltonian. Using the boson com- 
mutation relation for the cavity mode 



[A, A^] = AA^-A^A = 1 



one has 



— = -iu A - iKA f AA 
dt 



(12) 



(13) 



iK± / du>ai{uj) — iK2 I dua2(w) — i2n^ / du A** 0,2,(10) 



Using the boson commutation relations for the bath 
modes 



[ai(uj),aj(uj')] = S t j5(uj - w') 
[ai{w),aj{u)')\ = , 



(14) 
(15) 



one obtains the following equations for the bath modes 
a\{fjj), a 2 (uj), and a 3 (uj): 



da\(u>) 
dt 



-itoaiiuj) — in\A , 



(16) 
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IV. RESPONSE TO A CLASSICAL PUMP 



and 



da 2 (w) 
dt 



da 3 (uj) 
dt 



-iuja,2((jj) — ik\A 



(17) 



(18) 



Using the standard methods of Gardiner and Collett [15| , 
these equations yield the following equation for the cavity 
mode A driven by the incoming bath modes a": 

— = -iuj A - iKA^AA - 7 A - kA^AA (19) 
dt 

-i2 y /^e i ^^af{t) . 



where 



7 = 7i + 72 



(20) 



and the which in general can be complex, have been 
reexpressed in terms of the positive real constants ji and 
the phases <j>i according to 



Kx = 4/ — 



K 2 = \ — 



K3 




(21) 
(22) 
(23) 



Expressions for 72 and 73 in terms of linear and nonlin- 
ear resistance of the stripline (see Eq. (J2J) are given in 
Eqs. (|ljl4|) . (11315(1 . In addition, the methods of Gardiner 
and Collett [l5j yield the following relations between the 
outgoing bath modes a° uf , the incoming bath modes a\ n , 
and the cavity mode A 



af"(t) - a\ n (t) 



ar\t) - a 2 n (t) 



l A{t) , (24) 



>-A(t) , (25) 



Operated as a negative-resistance reflection amplifier, 
an intense sinusoidal field, called the pump, is delivered 
to the input port of the device. Signals having frequen- 
cies to either side of the pump, but lying within the band- 
width of the device, will be amplified. The linearization 
procedure is now carried out in which the signals entering 
the input port and the noise entering the loss ports are 
considered to be small compared to the pump. The first 
step is to calculate the classical response of the device to 
an intense pump in the absence of signal and noise. The 
solution is then used to calculate the linearized response 
of the device in the presence of signal and noise. 

In order to obtain the response of the device to a clas- 
sical pump in the absence of signal and noise one sets the 
incoming noise terms to zero 



a 2 " = 0, 

4" = o. 

The incoming pump is written as 

a in — 5«« e -i("pt+^i) 



(29) 
(30) 

(31) 



where b™ is a real constant, oj p is the pump frequency, 
and tpi is the pump phase. The outgoing field will also 
have an oscillatory time dependence of frequency lo p and 
can be written as 



(32) 



where b° ut may be a complex constant. Writing A as 



(33) 



where B is a positive real constant, the equations of mo- 
tion, Eqs. {EH and yield 

[i(w -w p )+7]S+^+ 73 )S 3 = -ly^b^e^+'t'z-^ 

(34) 

and 



b out _ b in _ i 0^_ Be -i(0 1 +^ s _V, 1 ) 



(35) 



Multiplying each side of the Eq. 1)34(1 by its complex con- 
jugate and introducing 



E = B 



(36) 



(26) 



In obtaining these equations a Markov approximation 
|l5j has been made such that the boson annihilation op- 
erators a\ n (t) satisfy the commutation relations 



{ar(t),a^(t')}=6 h3 6(t-t>) 
[ar(t),ant')]=0. 



(27) 
(28) 



one obtains 



E' 



2[(u> - u) p )K + 773] £, 2 , (^o-^p) 2 +7 



271 



K 2 



7! 



K 2 



7 3 







-E 

(37) 



This cubic equation will have one real solution and two 
complex solutions or three real solutions. For the case 
when two of the solutions are complex, the real solution 
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FIG. 2: The cavity-mode amplitude \B\, the reflection am- 
plitude | | , and the intermodulation gain Gi for van- 
ishing offset frequency uo — shown for sub-critical case 
b\ n = 0.56^, critical case tf™ = b^c, and above-critical 
case &i™ = 2&ic- In au cases K = — 10~ 4 u , 71 = O.Olwo, 
72 = I.I71, and 73 = O.OlK/y/3. For bf > V{q the response 
becomes multi- value function of frequency in some frequency 
range. 



is the physical solution. If there are three solutions, two 
will be stable and one unstable, and the device will ex- 
hibit bistability. Once E and, hence, B have been deter- 
mined from Eqs. i|37|) and i|36|) . the phase <\>b can be de- 
termined from Eq. I|34|) and the amplitude of the reflected 
pump can then be computed from Eq. 1(35}. In Fig. 2 
(a), (d), and (g) plots of B as a function of frequency 
for three different incoming pump amplitudes are shown. 
The frequency pulling of the cavity resonance is clearly 
seen in (d) and (g) as the incoming pump amplitude b™ 
is increased. Also plotted in Fig. 2 (b), (e), and (h) is 
the reflection coefficient |6°"*/6™| for the reflected pump 
amplitude as a function of frequency. If no power were 
absorbed by the cavity, the reflection coefficient would 
be unity. One sees a dip in the reflected power at the 
cavity resonance. As the incoming pump amplitude is 
increased, this absorption feature also shows frequency 
pulling, as can be seen in (e) and (h). 



A. Special operating points 

As a function of the pump frequency lu p , B will have 
the form of the distorted Lorentzian curve (see Fig . 2 
fa) , (d), and (g)) exhibited by Duffing oscillators pMErj 
|27| . The maximum of the response curve occurs when 
dE/dutp = 0. This condition yields 



uj - uj p + KE = 



(38) 



that is, the peak of the resonance curve is shifted by 
an amount KB 2 . The points of instability where the 



system will switch from one of the two bistable states to 
the other are located where doj p /dE — 0. This condition 
is satisfied when 

( 7 + 2 73 £) 2 - (K 2 + 7 2 )£ 2 - (cjo - oj v + 2KE) 2 . (39) 

When, in addition, d 2 oj p /dE 2 = 0, the two points of 
instability coalesce into a single point. The condition 
d 2 uj p /dE 2 = is satisfied when 



6(K 2 + jDE + 4[(wo - u P )K + 773 ] = 



(40) 



Large parametric gain is achieved at points where the 
slope of E with respect to oj p becomes infinite, but in 
order to remain stable it is desirable to operate the re- 
flection parametric amplifier near the critical point with 
parameters chosen so that the Duffing curve does not 
have a bistable region. It is a straightforward exercise to 
show that in order for the resonance curve to have a crit- 
ical point at which both Eqs. I|39|l and (|40(l are satisfied, 
one must have 



\K\ > V3 73 
At the critical point one has 

27 



and 



w 



E c 



K 



V3(\K\ - ^73) 



4 73 |if| + V3(X 2 + 7 3 2 ) 
K 2 - 3 7 f 



(41) 



(42) 



(43) 



The incoming pump amplitude required for operation at 
the critical point is given by 



(6? 



in\2 



7 3 (if 2 +7!) 



3^3 71(1^1-^73 



(44) 



Thus, the input power required for driving the system 
into the threshold of bistability (critical point) is in- 
creased in the presence of two-photon losses. More- 
over, when 73 exceeds the value of \K |/v3 the bistability 
regime becomes inaccessible (see Eq. I|41|l). 
When 73 = these reduce to 



E c = 



LOq — LO p 



(&T 



m\2 



2V37 

m\ 

4 



K 

W\ 

7 3 



3V3 7i|^l 



(45) 
(46) 
(47) 



In panel (d) of Fig. 2 the amplitude of the cavity mode 
as a function of frequency has been plotted for the case 
when the incoming pump amplitude is that of the critical 
pump amplitude. One sees that the line shape of the 
cavity mode is vertical at a point on the lower side of the 
resonance. The line shape of the reflected power is shown 
in panel (e) of Fig. 2. 
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V. LINEARIZATION 



VI. SOLVING THE LINEARIZED EQUATION 



A linearized analysis is now performed in which the 
incoming signal and the noise from the losses are regarded 
as small compared to the pump. To that end we write 



and 



<2j" = 5™g- i (^p t +'/'i) _|_ c » 



in ^ — iUpt 



„in in —iui v t 

a 2 — c 2 6 



in An^ — iu v t 



IjOut g—i(uj p t+tl>i) _j_ gOutg— iwpt 



out lout —iujpt i out —iujpt 

Ao — Co C ~ p Co C 



out uout„ — iujpt i „out „— iu)pt 

3 



(48) 



(49) 



(50) 



(51) 



(52) 



(53) 



(54) 



where J5, 6°"', 63"' and 63"* constitute the solution for the 
response of the system to a classical pump in the absence 
of signal and noise. The properties of this solution have 
already been discussed in Section 4. The c™, c™, c™, 
cf"*, Cj"*, Cg"', and a are regarded as small and will be 
kept only up to linear order. Substituting these into the 
equations of motion yields 



Introducing 

W = i(Lu -Lu p )+j + 2{iK + -/ 3 )B 2 , (59) 



2 -2i 



and 



V = (iK + l3 )B z e 



(60) 



(61) 



the linearized equation of motion can be written in the 
form 



da 
~dt 



Wa + Va) = F 



From this last equation one obtains 

^ + 2M{W)^ + (\W\ 2 -\V\ 2 )a = T(t) 



where 



Writing 



dP 

r(t) = — + W*F -VF\t) 



-XI 



(62) 



(63) 



(64) 



(65) 



the characteristic equation for the homogenous equation 
is given by 

A 2 -23f?HA + |W^| 2 - \V\ 2 = . (66) 
This has the two roots 

A - K(W) - y/W(W)- M 2 + \V\ 2 ) , (67) 



da 



■ — [i{ujQ — uop) + 7] a — 2(iK + -f 3 )B 2 a 

(iK + 73 ) J B 2 e- i2lfe a t - iy^e^ci 71 
iy/^e** 2 ^ - i2V^£e l ^ t+ ^ +< ^4™ , (55) 



or 



Ai = + y/% 2 (W)-\W\ 2 + \V\ 2 ) , (68) 



A = 7 + 273B 2 - yf (K 2 + 7 3 2 )S 4 - (coo - cj p + 2KB 2 ) 2 , 

(69) 



c out _ c m = _ iy /2^ e -i<f>i a } ( 56 ) 



c out _ c m = _j2 v ^ 3 "Be- l( ^ t+ ^ B+ ^ ) a . (58) 



Ai = ~f + 2~/ 3 B 2 + y (K 2 +7|)_B 4 — (co ~ u. p + 2KB 2 ) 2 . 

(70) 

The root Ao is zero when Eq. H39|l is satisfied, that is, one 
has critical slowing down at the points where the slope 
of E with respect to lo p is infinite. 
Introducing the Fourier transforms 



a(t) 



'2-k 



du>a(cj)e 



(71) 
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or 



d(f) 



Pa(*) = 



27T J_ 



27T 



C3 



1 f° 



r(t) 

Eq. S3 yields 

a(w) = 



27T J_ 



dwci (u>)e 



dwc3(a;)e' 



dwr(a;)e 



-iut 



(72) 



(73) 



(74) 



(75) 



w 2 - 2iwR(W) + (|W| 2 - |y| 2 ) 



(76) 



In terms of the roots of the characteristic equation this 
can be written as 



a(u>) 



I» 



(—icj + Xo)(—iiu + Ai) 



(77) 



where 



r(w) = -»v^n[(-iw + Wy^c^fw) - Ve-^cf^-w)] 

- ly^K-iuj + W*)e*4"(o;) - 7e-^ 2 4 nt (- w )] 

- Ve^* B+ *»^(u p -u)] . (78) 



A. The output field 
From Eq. (|24l) one obtains 

c° ut (w) = c™H - iy/^e-^a^u) . (79) 
Substituting Eq. I|77|l into this equation yields 



(-iu + X )(-iu + Ai)ci"(w) - iV27i e_i0ir (^) 



(—to; + Ao)(— iw + Ai) 



(80) 



{-iu + X )(-iu + Ai) - 2j 1 (-icj + W*) 



+ 



(—iuj + X )(—iu + Ai) 



4» 



(— + Ao)(— zcj + Ai) 
2V7i72(-^ + l^*) e - l (^-^) 
(— iu + X )(—iu + Ai) 

2V7iWe- i( * 1+ » 2) e<wt 
(-iw + A )(-^ + Ai) C2 WJ 

2 v / 27i7iS(-iw + H/*)e- 4 ^ 1 -^-^) inl 

c™(u p + u) 



( — ILU + Xo)( — lLU + Al) 

1+4>3+4>b) 



(—iu; + Xo)(—iuj + Ai) 



(Up - u) 



(81) 



The linearized solution to the equations of motion has 
now been obtained. We will now evaluate the properties 
of this solution for various kinds of inputs. 



VII. PARAMETRIC AND 
INTERMODULATION GAIN 

The parametric gain and the intermodulation gain are 
calculated by taking c™(w) to represent a classical signal 
at frequency u p + u. Setting all other signal and noise 
inputs to zero, Eq. 1|81[) yields the following power gain 
for the reflected signal: 



\ c out 



HI 5 



|cl«HI 2 

(-*w + Ao)(-*w + Ai)-2ti(- 



(c 2 + A 2 )(^ + A 2 ) 



(82) 



When this quantity becomes greater than unity one has 
parametric amplification of the signal. 

As seen from Eq. |(8TJ, a signal c in (—u) injected at fre- 
quency u p — u will generate an output signal at frequency 
Lo p + lo. This frequency conversion is quantified by the 
intermodulation conversion gain defined by 

_ ic^HI 2 

w = 



I4 n (-cj)| 2 
4 7l 2 M 2 

( W 2 + A 2 )( W 2 + A 2 ) 



(83) 



Since the output signal at uj p +lu is separated in frequency 
from the input signal, the measurement of the intermod- 
ulation gain is a particularly sensitive method for mea- 
suring the strength of the nonlinearities. We note that, 
even without power gain, devices capable of producing in- 
termodulation signals are useful as mixers. When u = 0, 
both the expression for Gs and the expression for Gi 
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will have AqA| in the denominator. As one approaches 
an operating point where the slope of E with respect to 
Up becomes infinite, both the parametric gain and the 
intermodulation conversion gain will diverge. Hence, it 
is near the instability points where the device can exhibit 
large gains. Panels (c), (f), and (i) of Fig. 2 show the 
behavior of the intermodulation gain as the pump am- 
plitude is increased from half critical (c) to critical (f) to 
twice critical (i) as a function of frequency. As depicted 
in panel (f ) at the critical point the intermodulation gain 
diverges. Above critical, as shown in panel (i), the inter- 
modulation gain diverges as one approaches the points of 
infinite slope on the resonance curve (g). 



VIII. NOISE SQUEEZING 

Because of intermodulation gain, a parametric ampli- 
fier can establish correlations [23] between the output at 
lu p + oj and lo p — uj. When delivered to a mixer whose 
local oscillator is phase-locked to the pump these correla- 
tions can result in noise fluctuations reduced below that 
which the mixer would see if the signal delivered to the 
parametric amplifier were, instead, directly delivered to 
the mixer. This noise reduction is called squeezing, and 
it can occur with either thermal or quantum noise 
We now obtain expressions that will allow one to calcu- 
late the degree of thermal or quantum noise squeezing. 



For such a calculation the c™, c 
are again treated as quantum mechanical operators sat- 
isfying commutation relations of the form Eqs. Q27J1 and 



-3 1 



The output of a mixer, operated in the homodyne 
mode in which the local oscillator frequency ujlo and 
the pump frequency lj p are equal and in which the input 
at the signal frequency ljlo + w and at the image fre- 
quency lolo ~ w are both regarded as signal, is given by 



4] 



out] 



(84) 



where 4>lo is the local oscillator phase. To evaluate the 
mean value and the power spectrum for the homodyne 
detector output it is necessary to specify the density ma- 
trix for the signal and noise entering the parametric am- 
plifier. Here we consider the case when these inputs con- 
sist of Nyquist noise. In this case one has 



(CM) = , 



(<ft(u;)c?V)> = ' 



—fiihojp 



Sij5(u - uj') , 



and 



Here 



(cf(u) C jV))=0 



(85) 
(86) 
(87) 



kBTi 



where is Boltzmann's constant and Ti is the absolute 
temperature of the bath for which c* n (w) is the incom- 
ing mode. We thus allow each of the baths to be at a 
different temperature. In writing Eq. I|86[) we have made 
the approximation that the frequencies lo of interest are 
small compared to u> p . 

Because Eq. I|81|l is linear in the c\ n (ui) it is evident 
that 



{I D (w)) =0 



(89) 



that is, the homodyne detector output consists of noise 
fluctuations with zero mean. Because of the boson com- 
mutation relations, one has 



(4(w)/2,(w')) = P(v)6(w - a/) 



(90) 



where P(u>) is the noise power spectrum of the homodyne 
detector output. Equation (|81|l can be rewritten as 

c? u 'M = Ai(w)cT(w) + B 1 ( W )cl" t (- W ) 
+ A 2 (u)c?{uj) + B 2 (uj)c^(-uj) 
+ A 3 {cj)c l 3 n {cjp + uj) + B 3 {Lj)cl n \ujp-uj) (91) 



where 



{-iuj + X )(-iu + Ai) - 271 (-iuj + W*) 
(— iu) + Ao)(— iu> + Ai) 



2 7 i Ve- a ^ 



(—iuj + \ a )(—iLU + Ai) 



(92) 



(93) 



A 2{u) = v . . , (94) 

(— iuj + \ q )(—ioj + AiJ 



B 2 (u) 



2V7i72^ e ~ i(01+<M 
(— iu + Aq)(— iu + Ai) 



(95) 



A 3{u) = — — . — — — , (96) 

(— iu + X )(-iuj + Ai) 



(-10J + \q)[ — 1(jJ + Ai) 

Substituting Eq. JSU into Eq. (JBU, evaluating 
(ily {u))Id (a/)) using Eq. JSSJl, and then reading off the 
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power spectrum using Eq. I|90[l . one obtains 



P(w) = \e- l4>LO A\{uj) + e l4,LO Bi{-u)\ 



+ |e^ L0 Ai(-cj) + e^ I ' B 1 *(w; 



+ |e^ A 3 (-w) + e~^ i0 B3*H| 



1 


R hi . 
_ g—piflUJp 




1 


1 


g M± p 




O — 02 Ti>bJ v 


\ 
1 


g — &2^p 




1 


2 T 


_ e -&swp 






2 r 


— g—PaUwj, 




1 


2 t 


_ g— foHup 



(98) 



This formula may be used to compute the noise-power 
spectrum for any local oscillator phase <Plo and any set 
of device parameters. It is useful to consider the case 
when there is no incoming pump field and the input field 
and loss baths are all at zero temperature. In this case 
the field reflected off the input port of the amplifier will 
consist of vacuum fluctuations, that is, 

c r'Mio) = o , 



and one obtains 



(99) 
(100) 
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FIG. 3: Examples of achievable degree of squeezing, P m i n (0) 
vs. bf/bTc- In all plots Ti = T 2 = T 3 = 0, K = 5lo , 
7i = O.OOOIojo, and the pump frequency uj v obtains its critical 
value, as given by Eq. 1431 The solid line represents the loss- 
less case, where 72 = 73 = 0. The dashed line represents the 
case of linear dissipation where 72 = 571 and 73 = 0, while 
the dotted line represents the case of nonlinear dissipation 
where 72 = and 73 = 0.5if/v3- 



A. Special cases 



P{w) = 1 



(101) 



This sets the vacuum noise level for the conventions we 
are using. As is shown in Fig. 3 and as will be illus- 
trated with specific examples in the next section, under 
suitable circumstances it is possible to obtain reflected 
signals whose noise-power spectrum P(uj), for certain lo- 
cal oscillator phase settings, is less than 1. Such signals 
are said to be squeezed below the vacuum noise level. 
Figure 3 shows the minimum value of P(lo) as a function 
of the amplitude of the incoming pump 6™ when the 
bath temperatures are all zero, the strength of the Kerr 
nonlinearity is chosen to be K = 5loq, and the coupling 
strength of the signal to the cavity mode is taken to be 
71 = 0.0001^0- The pump frequency has been chosen to 
be that of the critical pump frequency. The solid line is 
the lossless case when 72 = 73 = 0. In this case complete 
squeezing P(0) — is possible when the incoming pump 
power is at the critical value, that is, = 1. The 

dashed curve represents the case when the linear dissipa- 
tion is given by 72 = 571 and the nonlinear dissipation 
73 is zero. In general, the presence of linear dissipation 
reduces the amount of achievable squeezing. The dotted 
line depicts the case when the linear dissipation is zero 
and the nonlinear dissipation is given by 73 = 0.5K/V3. 
Nonlinear dissipation also tends to reduce the amount of 
squeezing that can be produced. However, nonlinear dis- 
sipation can produce some squeezing in the absence of a 
Kerr medium, as we discuss below. 



It is instructive to evaluate Eq. (|98|l for a few specific 
cases. In particular, the power spectrum at uj — reduces 
to 



P(0) 



>AI(0) 



'Bi(0)| 2 coth 



i4>i>oA*(0) + e i<AiO B 2 (0)| 2 coth 
i(f,LO -B 3 (0)| 2 coth 



/ fkUp 
\2k B T! 



2k B T 2 



2k B T 3 



(102) 



Further simplification results from considering the case 
when the internal losses are all zero. In this case 72 = 
and 73 = 0. Thus, the terms in Eq. 1102fl corresponding 
to internally generated noise vanish and one obtains 

P(0) = \e- i<t,L °A{(0) + e l ^° J B 1 (0)| 2 coth ( 

\2k B Ti 



Writing 



one has 



A 1 {Q) = \A l {Q)\e i ** , 
Si(0) = |Pi(0)|e^ fl , 



P(0) = ||Ax(0)| + l^^le^fcoth (^g^ 



(103) 

(104) 
(105) 

(106) 
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where 



IX. CONCLUSIONS 



V> = 24>lo - 4>a + 4>i 



(107) 



P(0) is maximized when the local oscillator phase 4>lo is 
chosen so that 



1 



In this case one obtains 



P max (0) = (1^(0)1 + |Bi(0)|) 2 coth 



/ hujp 
\2k B T x 



(108) 
(109) 



P(0) is minimized when the local oscillator phase §lo is 
chosen so that 



i-0 



1 . 



In this case one obtains 



P™„(0) = (|Ai(0)|-|Bi(0)|) a coth 



ftWp 

2fc B Ti 



(110) 
(111) 



From Eqs. (fTUTjl . l(TU5|) . and lfTTU|) it follows that the lo- 
cal oscillator phase at which the spectral density P(0) is 
minimized differs by 7r/2 from the local oscillator phase 
at which the spectral density is maximized, that is, the 
signal components minimizing and maximizing the spec- 
tral density are in quadrature. It is straightforward to 
show that 

|^-x(0)j 2 - |^i(0)j 2 = 1 . (112) 
From this equation and Eqs. I|109|) and (|lllfl . one obtains 



P m ax(0)Pmin(0) = COth' 



( hujp 



\2k B T 1 



(113) 



Hence, in the case of no loss, the degree of amplification 
and the degree of deamplification of the noise are the 
same. For this case one also has 

2^KB 2 



\Bi(0)\ 



(«0 -uj p + 2KB 2 ) 2 + 7 f 

Note that |i?i(0)| can be made arbitrarily large by choos- 
ing the pump frequency lo p and the pump amplitude B 
such that the denominator goes to zero. From this it fol- 
lows that P m ax can be made arbitrarily large and P m in 
can be made arbitrarily small. In practice, higher order 
terms responsible for pump depletion will limit how big 



K 2 B* 



(114) 



We have presented an analysis of a cavity parametric 
amplifier employing a Kerr nonlinearity but which also 
possesses a two-photon loss. We have obtained expres- 
sions for the pump amplitude inside the cavity and the 
reflected pump amplitude for the case when pump satura- 
tion can be neglected. We have obtained expressions for 
the classical gain and the intermodulation gain. These 
expressions are useful for determining model parameters 
from experimental data. We have found that in the pres- 
ence of two photon losses the injected power required for 
driving the system into the bistable regime is increased. 
Moreover, this regime becomes inaccessible when 73 ex- 
ceeds the value of |i<r|/\/3. 

We have also obtained expressions from which one can 
compute the degree of squeezing that the device exhibits. 
Both the linear and nonlinear loss tend to degrade the 
amount of squeezing that can be achieved, although even 
without a Kerr nonlinearity a modest amount of squeez- 
ing can be achieved by the two-photon loss. 



APPENDIX A: NONLINEAR KINETIC 
INDUCTANCE IN TRANSMISSION LINE 
RESONATOR 

Consider a lossless linear transmission line with length 
I extending along the x-axis. Let q (x, t) be the charge 
density per unit length and define |l9j | 



(x, t) = / dx'q (x' , t) . 



(Al) 



Thus q = —dQ/dx, and the voltage across the trans- 
mission line is given by 



V(x,t) 



i_dQ 

C dx'' 



(A2) 



where C is the capacitance per unit length along the 
transmission line, whereas the current is given by 



Noise from the losses and gain saturation will limit 
the maximum degree of noise squeezing to be below that 
calculated here. This is demonstrated in Figure 3, where 
the achievable degree of squeezing P m i n (0) is plotted as a 
function of for some examples in which the linear 

or nonlinear loss is taken to be nonzero. 

The presence of a two-photon loss allows for squeez- 
ing even in the absence of a Kerr nonlinearity [23j . Set- 
ting K — 0, 72 = 0, and ujq — uj p , the greatest de- 
gree of squeezing of the output field of the cavity occurs 
when 7! = 373-B 2 . At this operating point P(0) = 2/3 
and the local oscillator phase must be adjusted so that 
cos(20 lo -20b-20i) = 1. 



I{x,t) 



dQ 
dt ' 



(A3) 



The Lagrangian C of the system reads 



1 

C=^Jdx [LI 2 - CV 2 ] 



dx 



dQ 
dt 



(A4) 



I (09. 

C [ dx 
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where L is the inductance per unit length along the 
transmission line. The open ends at x = and x = I 
impose boundary conditions of vanishing current. 

We assume the case of a nonuniform transmission line, 
where both C and L may depend on x. Moreover, the 
inductance L depends on the current / according to Eq. 
as a result of nonlinear kinetic inductance. 

As a basis for expanding Q (x,t) as 



Q ( x > *) = E qn (*) Un ( x ) ' 

71 

we use the solutions of 



(A5) 



Tx {hit) = (A6) 

with the boundary conditions of vanishing current 



U n (0) = U n (I) = 0. 



(A7) 



We assume that the functions u n (x) are chosen to be 
real. 

For the case of a uniform transmission line resonator, 
where both C and L are independent on x, one has 
u) n = mr/l\/LC, where n is an integer. In the nonlinear 
regime, however, such an equally spaced spectrum may 
lead to strong inter-mode coupling, where harmonics and 
sub-harmonics of a driven mode excite other modes. In 
the present work we assume that such inter-modes ef- 
fects are avoided by employing a nonuniform resonator. 
This allows us to consider only the mode in the resonator 
which is driven externally. 

Using expansion (|A5|) 



n m q 

~ \ EE q,iqm J dx 



(A8) 



1 du n du m 



C dx dx 



AC, 



where 



A£ = in' in" in" 'in" 



TV ,n" ,n' ,n" 



x J dx u n 'U n "U n 'f'U n 'ff AL. 





(A9) 



We first treat the linear part. Consider equation (|A6|I 
for u n multiplied by u m and equation \A&\ for u m mul- 
tiplied by u n 



d ( 1 du n \ > 
Um dx \C~dx~ ' = ^^nLoUnUm, (A10) 



1 du m 
dx \C dx 



-uj m Lou n u v 



(All) 



Subtracting 



1 du n 



1 du, 



dx \ Um C dx U " C dx 



and integrating from x — to x — I one obtains 



(w, 2 „ - Wn) / dxL u n u m = 0. 



(A12) 



(A13) 



In general, it can be easily shown that the spectrum of 
a finite one-dimensional resonator having vanishing cur- 
rent boundary conditions is non-degenerate. Thus, by 
requiring that the functions u n (x) are normalized one 
obtains 



dxL u n u m = 5 n 



(A14) 



Moreover, integrating by parts and using (| A6|) and the 
boundary conditions one obtains 



1 1 



C dx dx 



u) n / dxL u n u m = Lu n S nm . (A15) 











Thus 



C = \ E Un - ^nln) + A£. 



The Euler-Lagrange equation is given by 



DC 



dt \ dq n J dq 



dC 



(A16) 



thus 



9n+ ^„ + -^^0. (A17) 
The variable canonically conjugate to q n is 



_ dC_ _ . dAC 

Pn — qn T" ~ 

oq n oq n 
The Hamiltonian is given by 



(A18) 
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H = ^2p n q n - C 

n 



(A19) 



a fdACV 2 2 



AC. 



To first order in AC (or in AL) 



where 



H = H a + V, 



n = \^2{p 2 n+ w n«n) » 



(A20) 



(A21) 



and 



2/ c 2 



Pn'Pn"Pn ,/ 'Pn'' 



(A22) 



x dx u n iu n "U n "'U n "// AL. 



To quantize the problem the variables g„ and p n are re- 
garded as operators satisfying the following commutation 
relations 



[qn,Pm] = QnPm - P m q n = ifrSn 



[q n ,q m ] = [Pn,P m ] = 0. 



(A23) 



(A24) 



The Boson annihilation and creation operators are de- 
fined as 



An = —= ^/uJ^q n H —Pn 

V2h \ 



At. 



-iuj n t 



2h 



(A25) 



(A26) 



The inverse transformation is given by 



(A27) 



The commutation relations for the operators A n and 
A n are derived directly from l|A23() and (|A24|1 



[A n , A m ] — 8 n ,; 



[An, Am] — [A n ,A m ] — 0. 



(A29) 



(A30) 



Using (|A27|) and (|A28|) . the Hamiltonian l|A21() can be 

expressed as 



Ha = Y J r ^ J n [A n An + \ 



(A31) 



The current operator is given by 



n 

The voltage operator is given by 



A n e'^ 1 ) u n (x) 
(A32) 



du r 
dx 
(A33) 

The terms p n 'Pn"Pn'"Pn"" contain in general terms os- 
cillating rapidly at frequencies on the order of the fre- 
quencies in the resonator spectrum. In the rotating 
wave approximation (RWA) these terms are neglected 
since their effect on the dynamics on a time scale much 
longer compared to typical oscillation period is negligi- 
bly small and only stationary terms remain. Thus in the 
expression of V only terms of the type p n >P n " contain 
stationary terms, which are given by 



2 2 
Pn'Pn" 



nl(l + 2AlA n ,) (l + 2Al,A n 



2 2 

(A34) 

The constant term can be disregarded since it only 
gives rise to a constant phase factor. Moreover, the terms 
A n ,A n ' and A^ n ,,A n » that give rise to frequency shift can 
be absorbed into TCq. Thus in the RWA the perturbation 
V contain only terms of the type A n ,A n ' A n „A n " 



V = ^J2 K "'{ A l A "') 2+h E K'n»Al,An,Al,A n », 
n' n'^n" 

(A35) 

where 



fllOr, 



l >, l = lx \-^(Aie^ t -A n e-^) 



(A28) 



hup- 

K n , = — jf I dx u%AL, (A36) 
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and 



3tkJ T , 



P 



dx U„tU„r,AL. 



(A37) 



APPENDIX B: NONLINEAR LOSSES 
ASSOCIATED WITH THE KINETIC 
INDUCTANCE 

Here we derive expressions or the linear and nonlinear 
loss coefficients 72 and 73 in terms of the parameters 
that characterize the resistive loss associated with the 
kinetic inductance. This is accomplished by obtaining 
an expression for the rate of energy loss in the cavity in 
terms of the model parameters and comparing it with the 
expression for the power dissipated due to the resistance 
associated with the kinetic inductance. 

The equation of motion for the resonator Hamiltonian 
H r (Eq. ©) 



^— = [H r ,H] 



(Bl) 



yields 

dH r 



dt 



= —ihujQ / du)[KiA^ ai(w) — n*a\(uj)A] 



- ifruio J du>[K 2 A^ a 2 (ui) — K^a^^A] 

- 2ihuj J dw[n 3 A^ A^ a 3 {w) - K%a\{w)AA] 

- 2ihK J dw[K 3 A^ A^ a 3 (u>) - K%a\{w)AA] 

- ihK J du)[K X A^ A^ Aa x {u) - K\a\{u)A^ AA] 

- ihK J duj[K 2 A^ A^ Aa 2 {io) ~ n* 2 a\{ij)A^ AA] 

- 2ihK J dcj[K 3 A ] A ] A j Aa 3 (uj) - Kla\(u)rf AAA] 

(B2) 

As was shown by Gardiner and Collett 0, 0] , the equa- 
tions of motion Eq. Ijl6(l through ljl8|l for the baths can 
be integrated to yield 



1 



V2tt 
Similarly 

1 



dio ai(w) = a\ n (t) - i\j-K*xA{t) . 



/2tt 



duo a 2 {ijj) = a 2 n (t) - i\j-n 2 A(t) 



and 



duj a 3 (to) = a 3 n {t) - i\j -n%A{t)A{t) 



(B3) 



(B4) 



(B5) 



These expressions can be used to eliminate the a n (uj) 
from Eq. 1B2(I . Evaluating the expectation value of 
Eq. 1B2|I with respect to a state in which all the bath 
modes are in a vacuum state, that is, 



<(t)|0>=0, 



(B6) 



one obtains 



dH r 
~~dt 



= -2M,[7i<^) + TS<AU) + 273<A+AUA)] 

- 2hK[j 1 {A^A^AA) + l2 {A^A^AA) 

+ 2r/ 3 (A^A^AA) + 273^ A* A* A A A)] 

(B7) 



If luq ^> K then as long as the mean-field is not too large 
one has, to a good approximation, 



dH r 
dt 



\ = -2^ [ 7l (AU) + 72 (AU) + 273(AUUA)] 

(B8) 



This is an expression for the rate with which energy is 
lost from the cavity. 

The power dissipated within the cavity is given by 



P = / dx (: RI 2 



(B9) 



where the ": :" denotes normal ordering. Using Eq. J5J, 
this can be written as 

P= [ dx R (: I 2 :) + / dx Ai ?( : /4 : ) ■ ( B1 0) 

Jo 4 Jo 

where we have allowed for the possibility that i? an d Ai? 
may be functions of the distance x along the resonator, 
as would be the expected case if the composition and 
shape of the transmission line cross section varies with x. 
Substituting Eq. (|A32|I into this equation yields 

P = hu n [ dx R ul(AlA n ) (Bll) 
Jo 

In order to make the comparison of this equation with 
that of Eq. I|B8[) we consider the unloaded cavity case 
when the coupling through the signal port of the cavity 
is set to zero, that is, 



71 =0 



Setting 



P = - 



dH r 
~dt 



(B12) 



(B13) 
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and keeping in mind that presently wo and u n both de- 
note the resonance frequency of the same selected mode, 
one has 



and 



73 



72 = r 



dx uiRo 



(B14) 



dx uiAR . 



(B15) 
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